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ABSTRACT 

We investigate the co-evolution of black-hole-accretion-rate (BHAR) and star-formation-rate (SFR) 
inl.5<2:<2.5 galaxies displaying a greater diversity of star-forming properties compared to previous 
studies. We combine X-ray stacking and far-IR photometry of stellar mass-limited samples of normal 
star-forming, starburst and quiescent/quenched galaxies in the COSMOS field. We corroborate the 
existence of a strong correlation between BHAR (i.e. the X-ray luminosity. Lx), and stellar mass (M*) 
for normal star-forming galaxies, although find a steeper relation than previously reported. We find 
that starbursts show a factor of 3 enhancement in BHAR compared to normal SF galaxies (against a 
factor of 6 excess in SFR), while quiescents show a deficit of a factor x5.5 at a given mass. One possible 
interpretation of this is that the starburst phase does not coincide with cosmologically relevant BH 
growth, or that starburst-inducing mergers are more efficient at boosting SFR than BHAR. Contrary 
to studies based on smaller samples, we find the BHAR/SFR ratio of main sequence (MS) galaxies 
is not mass invariant, but scales weakly as ]VfO-43±o.09^ implying faster BH growth in more massive 
galaxies at z 2. Furthermore, BHAR/SFR during the starburst is a factor of 2 lower than in MS 
galaxies, at odds with the predictions of hydrodynamical simulations of merger galaxies that foresee 
a sudden enhancement of Lx/SFR during the merger. Finally, we estimate that the bulk of the 
accretion density of the Universe at z ~ 2 is associated with normal star-forming systems, with only 
^ 6(±1)% and ~ 11(±1)% associated with starburst and quiescent galaxies, respectively. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: nuclei — galaxies: starburst 


1. INTRODUCTION 

The self-regulation of host galaxies and their central 
supermassive black holes (BH) is often invoked to ex¬ 
plain the local scaling relations between BH mass and 
galaxy properties, (such as the stellar mass or the ve¬ 
locity dispersion of the bulge; Magorrian et al. 1998, 
Ferrarese & Merritt 2000; Aller & Richstone 2007). As a 
consequence, determining the relative roles of the various 
processes that drive the co-evolution of BHs and galax¬ 
ies has emerged as a key goal of current astrophysics 
research. 

While the role of violent mergers is still considered as 
an important triggering mechanism for BH growth (ob¬ 
served as active galactic nuclei; AGN), in particular in 
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very luminous quasars (Hopkins et al. 2008, Treister et 
al. 2012), the emerging observational framework sup¬ 
ports the idea that secular processes are responsible for 
both star formation and SMBH growth in a large ma¬ 
jority of galaxies displaying moderate nuclear activity 
(Cisternas et al. 2011, Kocevski et al. 2012). As such, 
it is now clear that processes other than galaxy merg¬ 
ers, such as gas turbulence or disk instabilities, play an 
important role in triggering BH accretion (Bournaud et 
al. 2012). However, with many studies finding no clear 
correlation between instantaneous BH growth and star 
formation (SF; Silverman et al. 2009, Shao et al. 2010, 
Rosario et al. 2012, Mullaney et al. 2012a, Azadi et al. 
2014), actually measuring the contribution of these dif¬ 
ferent triggering processes to the BH growth budget has 
proven a significant challenge. 

A key step forward in this area has been achieved by 
studies that have adopted X-ray stacking analyses to de¬ 
rive the average BH growth rate of large samples of galax¬ 
ies (i.e., including undetected, low accretion rate BHs). 
In doing so, these studies have demonstrated that the 
average BHAR is tightly correlated with both M* and 
SFR (Mullaney et al. 2012b, hereafter M12; Ghen et al. 
2013) and mimics the SFR-M* relation of normal star¬ 
forming galaxies (Elbaz et al. 2007, Daddi et al. 2007, 
Rodighiero et al. 2014). However, these previous studies 
did not discriminate between different types of galaxies, 
whereas it is known that different levels of SFR are trig¬ 
gered by different events which may also trigger different 
levels of BH growth. Moreover, these studies did not 
consider BHAR in quiescent galaxies, leaving their con¬ 
tribution unexplored. 
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To constrain the relative BHAR in galaxies display¬ 
ing different levels of star-forming properties (and thus, 
in different evolutionary stages), we employ X-ray stack¬ 
ing to investigate how the average X-ray luminosity (a 
proxy for BH growth rate) varies across bins of stel¬ 
lar mass and specific SFR (sSFR) using large statisti¬ 
cal stellar mass-selected samples. Exploiting the unique 
depth/area combination of the COSMOS field (Scoville 
et al. 2007) means we avoid the small number statistics 
that prevented previous X-ray stacking studies to separ¬ 
ate their samples in terms of sSFR. Here, we focus on 
the 1.5 < z < 2.5 epoch when most of today’s stars were 
formed; see also Delvecchio et al. (2015, subm.). 

Throughout the paper we use a Chabrier (2003) initial 
mass function and = 71 km s“^, Ha = 0.73, Hm = 
0.27. 

2. SPLITTING THE SFR-M* PLANE AT 1.5 < z < 2.5 

2.1. Normal star-forming galaxies 

We use the K-band selected sample of 1.5 < z < 2.5 
star-forming galaxies down to Kg as < 23 in the COS¬ 
MOS field (McCracken et al. 2010) selected according to 
the BzK criterion designed to select star-forming galaxies 
at these redshifts (the star-forming BzK, hereafter sBzK; 
Daddi et al. 2004). Stellar masses have been computed 
following the same procedure as in Daddi et al. (2007). 
We have considered the potential issue of AGN contami¬ 
nation in the computation of stellar masses (in particular 
in the highest mass bin): we verified that the AGN stel¬ 
lar masses computed with a spectral decomposition that 
accounts for a dusty torus emission, provides on average 
consistent results with M* classically derived by a pure 
stellar component (no offset and a dispersion of ~0.25 
dex for both type 1 and type 2 AGN).SFRs of these sBzK 
galaxies are estimated from the UV rest-frame luminos¬ 
ity corrected for dust extinction, with reddening being 
inferred from the slope of the UV continuum as in Daddi 
et al. (2007). However, to factor-in the potential con¬ 
tamination from the type 1 AGN to the UV emission, we 
will use the average SFR(IR-I-UV) derived from stacking 
the sBzK sources on the PACS /Herschel maps in differ¬ 
ent mass bins (as in Rodighiero et al. 2014). 

Limiting our analysis to M* >1O^°M0, where the sam¬ 
ple is complete (see Ilbert et al. 2010), results in a final 
sBzK sample of 16,435 galaxies. See black points in Fig¬ 
ure 1, top panel. 

2.2. Starburst sources 

For our starburst sample, we start from the Gruppioni 
et al. (2013) sample selected from PACS/Herschel PEP 
observations in the GOSMOS field (Lutz et al. 2011) 
at lOO/xm and 160/xm.We refer to Berta et al. (2010) 
and Gruppioni et al. (2013) for a description of the data 
catalogues. Stellar masses have been derived by fitting 
the broad-band SEDs of our sources using a modified 
version of MAGPHYS (Da Cunha et al. 2008, Berta et 
al. 2013). Such M* are consistent with those derived 
with the formalism of Daddi et al. (2007, zero offset and 
a dispersion of ~0.3dex). The SFRs of these galaxies are 
derived from the total IR luminosities. 

To classify a sample of starburst galaxies lying above 
the MS, we have applied the same criterion as Rodighiero 
et al. (2011, i.e., sSFR 0.6 dex above the MS) over the 
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Fig. 1.— Top panel: The mass-SFR relation. We re¬ 
port the sBzK sample with their SFR(UV) (black points), 
while the corresponding average SFR(IR-hUV) derived from 
the PACS stacking analysis are reported as red filled circles 
(Rodighiero et al. 2014). The linear fit provides an MS re¬ 
lation as log[S'Fi2]=0.89(di0.09)xlog[M*]+14.02. We also show 
the starburst sample (blue points) with the corresponding aver¬ 
age SFR(IR-I-UV) in the two mass bins (errors on average SFR 
are the standard deviation of the SFR distribution within each 
mass bin). Middle panel: stellar mass—average Lx relation at 
1.5 < 2 < 2.5 from our X-ray stacking analysis. Shaded colour re¬ 
gions are the results for the three samples: starburst (blue), main 
sequence (red) and quiescent sources (green). We provide a fit 
for the MS sample (solid red line, log[Lx]=l-40(ib0.09)xlog[M+]- 
14.11). The open cyan squares are the results of M12 in GOODS- 
S, scaled to the same K-correction adopted in this work. The 
dashed and dotted black lines represent the average Lx of de¬ 
tected and stacked X-ray sources in the MS, respectively. Bot¬ 
tom panel: Lx/SFR relation for MS (red) and starburst (blue). 
The dotted line is a linear fit to the relation for the MS sam¬ 
ple (log[Lx/SFR]=0.43(di0.09)xlog[M+]-|-35.6) and the open cyan 
squares are again the results of Mi2. In each panel the width of 
the stacked data along the x-axis represents the RMS of the mass 
distribution in each bin. 

whole 1.5 < z < 2.5 range. This selection provides a total 
number of 231 candidate starbursts with log(M*)> 10.27 
Mq. See blue points in Figure 1, top panel. 

2.3. Quiescent galaxies 

To complete the study of galaxies populating the SFR- 
M* plane, we include the complementary sample of pas¬ 
sive BzK (pBzK). A description of the pBzK sample (in- 
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eluding photo-z) in the COSMOS field is presented in 
Strazzullo et al. (A&A submitted). M* have been de¬ 
rived as for sBzK parent sample and SFR derived via 
SED fitting are on average consistent with zero. How¬ 
ever, as well as purely quenched, red and dead galaxies, 
the pBzK sample also consists of galaxies with low SFR. 
The adopted selection provides a total of 2559 quiescent 
sources with log(M*)> 10.4 Mq. 

3. X-RAY DATA 

The central 0.9 deg^ of the COSMOS field has been 
imaged by Chandra for a total exposure time of 1.8 Ms, 
down to a depth of ^ 2 x 10“^® erg cm“^ s“^ in the 
soft band (Puccetti et al. 2009). A main catalog of 
1761 X-ray point sources detected in at least one band 
(0.5-2, 2-8, or 0.5-8 keV) has been published by Elvis et 
al. (2009), together with optical identifications, optical 
magnitudes and spectroscopic and photometric redshifts 
from different spectroscopic campaigns and/or from SED 
fitting estimates (Civano et al. 2012). 

We cross-correlated the different samples described 
above with the Civano et al. (2012) counterparts cata¬ 
log (using a matching radius of 1”). We considered only 
objects falling in the Chandra area (~ 50% of the total 
samples), recovering a total of 317, 19 and 47 matches for 
the sBzK, PEP and pBzK samples, respectively. The C- 
COSMOS exposure can be considered complete down to 
L(0.5-8 keV)=3xl043 gj.g g-i in the entire redshift range 
explored in this work. As such, these X-ray detections 
probe the moderate to high luminosity tail of the AGN 
population {Lx> 10^^ ^“"^^ erg s“^), complementing the 
deeper exposures in the GOODS fields (see M12). 

3.1. Stacked and average fluxes 

In order to study the average X-ray properties of 
sources not individually detected at the limit of the C- 
COSMOS image, we used the CSTACK tooO (v3.1, 
T. Miyaji). To reduce to a minimum the contamina¬ 
tion of detected sources in the stacked signal, we consid¬ 
ered those galaxies without a known X-ray source within 
10” from the centroid. This led to the exclusion in the 
analysis of about 5-10% of sources in the various mass 
bins. The total sample of stacked objects overlapping 
the Chandra area consists of 6393, 73 and 1008 objects 
in the sBzK, PEP and pBzK samples, respectively, at 
M* > IO^^Mq. The numbers of galaxies in each mass 
bin are reported in Table [TJ The count rates in the 0.5-2 
(soft, S) and 2-7 keV (hard, H) bands are computed by 
averaging the background subtracted count-rates mea¬ 
sured at each input position in the stacked sample (see 
Cimatti et al. 2013). The CSTAGK tool provides the 
results of a bootstrap re-sampling analysis, which uses a 
large number (500) of re-sampled stacked count rate. We 
take the median of the distribution as the best estimate 
of the stacked count rate, and assume the 5th to 95th 
percentile range as the associated confidence interval. 

From the ratio of the soft and hard bands we derived 
the Hardness Ratio, HR=(H-S)/(H-|-S). The values of 
HR are the average of the distribution of the values de¬ 
rived by resampling (thousands of times) the distribution 
of the corresponding stacked count rate.The HR for the 
stacked sBzK sample as a function of mass are shown 

http://cstack.ucsd.edu/cstack/ 
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Fig. 2.— Top panel: Stellar mass - HR relation at 1.5 < z < 2.5 
for the MS-sBzK stacked sample (filled black squares). The dotted 
line is a linear fit to the distribution {HR=0.31xlog[M*]-3.3). The 
vertical red dashed lines show the marginal constraints that we can 
derive for the starbursts, basically consistent with the MS. The 
values for the detected sources are also reported (green squares). 
Middle panel: Derived values for two different assumptions of 
the X-ray spectral slope F. Bottom panel: Final K-corrections for 
the two different values of F turn to be almost constant with mass. 
We adopt a constant value of r=1.4 and N/f=22. 

in the upper panel of Figure 2. Although a trend of in¬ 
creasing HR with mass can be seen, the relation is also 
consistent with being flat. We therefore applied the same 
conversion factors from a single power-law in all mass 
bins to derive fluxes and luminosities. An average value 
of HR^O (hence Nh ^ 23) is similar to results of stacked 
studies of much smaller samples of undetected galaxies 
in the M*-SFR plane (e.g. Olsen et al. 2013 reported an 
HR^0.09 for their SF sample). 

When combining the detected and the stacked X-ray 
signals to calculate the total average X-ray flux (see mid¬ 
dle panel of Fig.l) we first resample the corresponding 
distributions of flux values (assuming a Gaussian error 
distribution for the detected fluxes) and then derive the 
distribution of averaged X-ray fluxes using: 

Tz -I- ristacked * ^stacked 

^ ave — ~ 7 \^) 

^detected “r ^stacked 

where ridetected and ristacked are the numbers of X-ray de¬ 
tected sources and of stacked sources, respectively, while 
detected (with i = l.ridetected) and ^stacked are the cor¬ 
responding X-ray fluxes. As above, we take as our best 
estimate of the combined stacked-Fdetected X-ray flux 
the median of the distribution of the Eave values with 
a 95% confidence interval derived from the same distri¬ 
bution. When the stacked signal alone is not significant 
(i.e. consistent with 0 within 2a, as for the starburst 
and quiescent samples) we assume as a measure of the 
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lower and upper limits of the stacked+detected X-ray 
fluxes (or luminosities) the corresponding 5th and 95th 
percentiles. The main reason for resorting to stacking in 
our work is to estimate integrated SFR and Lx outputs 
over well defined galaxy population, rather than typical 
average quantities for each object in the sample. 

Absorption-corrected rest-frame 2-10 keV fluxes (and 
corresponding luminosities) are derived from the count 
rates in the observed hard 2-7 keV band by assuming 
an absorbed power-law with r=1.4 with the relevant 
conversion factors (as derived by WEBPIMMS) and K- 
corrections. A different choice of P (e.g. P = 1.8) would 
provides a higher normalisation for Lx, but which is 
still consistent with the results obtained for P = 1.4 
and therefore would not change our main results. Sim¬ 
ilarly, applying a mass-dependent K-correction derived 
from the HR with the same power-law and varying col¬ 
umn density would only change Lx by a maximum of 
15% (lower panels of Figure 2). Table [T] lists the 2-10 
keV X-ray luminosities for all the stacked samples. 

We have verified that the average contribution to the 
final (stacked-bdetected) Lx from star-formation is al¬ 
ways lower than 10%, by using the SFR-Lx relation of 
Vattakunnel et al. (2012), i.e. within the uncertain¬ 
ties of our data, and thus it does not affect our results. 
We note that the fractional contribution of SFR to the 
stacked signal is much larger (~50%), but it does not 
affect significantly Pave in eq. 1. 

4. RESULTS AND DISCUSSION 

The average Lx (stack + detections) vs. M* for the 
MS, starburst and quiescent samples is shown in the mid¬ 
dle panel of Figure 1 (with the relation derived separately 
for X-ray detected and stacked sources). Our analysis 
shows a robust Lx — M* correlation (the hidden AGN 
MS of M12) for normal star-forming galaxies, confirm¬ 
ing its existence also at higher M* This is a superlin- 
ear relation with a slope of 1.40(±0.09) in logarithmic 
units. The presence of this sequence, together with the 
SFR-M* MS, has been interpreted as a footprint of the 
same secular processes that feed and sustain both the 
SFR and the BHAR from the galaxies gas reservoirs. 
We note, however, that the relation observed in this 
work is steeper than that reported by M12 (who found 
a slope of 1.05±0.36, consistent with their larger errors) 
at the same redshift. In particular, our increased statis¬ 
tics at higher masses provides robust constraints on this 
slope. Interestingly, starburst galaxies show higher levels 
of BHAR, with respect to MS galaxies at fixed M*. By 
contrast, quiescent galaxies have a deficit of an average 
factor 5.5 in terms of Lx with respect to the MS, sug¬ 
gesting that they belong to an evolutionary phase where 
the gas fuelling the central engine (and the SFR) is al¬ 
most exhausted (or the accretion is simply temporarily 
suppressed). Given the number densities and average Lx 
of starburst, normal star-forming and quiescent galaxies 
in the mass bin 10.75 < log{M^,) < 11.25 (where all sam¬ 
ples are complete, spanning the X-ray luminosity range 
41.8 < log{Lx) < 43.6 ), we argue that the accretion ac¬ 
tivity during the starburst phase is just ^ 6(±1)% of the 
cosmic integral BHAR at this redshift, while that of qui¬ 
escent sources is ^ 11(±1)%. The bulk of the accretion 
density (~ 83(±1)%) of the Universe is then associated 
with normal star-forming systems (similarly to what is 



Fig. 3. — Comparison of our results with predictions of hydro- 
dynamical models that simultaneously follow star formation and 
the growth of black holes during a galaxy—galaxy merger (Di Mat- 
teo et al. 2005, Hopkins et al. 2012). a) To define the simulated 
merger/starburst phase we report the SFR history of the models 
normalized to the pre-merger phase. The temporal window where 
(SFR)/(SFR pre-merger) is larger than a factor 4 for all models is 
marked by the two vertical dot-dashed lines. b) Evolution of the 
X-ray luminosity with time, c) Evolution of the X-ray luminosity 
to SFR ratio, as a function of time. Two different values of the 
virial velocity of the collisional galaxies are considered in the Di 
Matteo et al. (2005) model (?j=160km/s black lines, ?j=320km/s 
magenta lines). Our results are presented as a red rectangle for the 
pre-merger phase (corresponding to our MS normal star-forming 
sources) and as a blue region for the merger phase (i.e. our star- 
burst sample). The width of Lx /SFR from our data spans the full 
range covered in all mass bins (see Fig. ^ bottom panel). 

found at z < 1.2, Georgakakis et al. 2014). 

During the starburst phase, where SFR is enhanced by 
a factor of on average (Fig. 1, top panel), the average 
Lx is also enhanced, but by a smaller amount (a factor 
of ^3, Fig. 1 middle panel). The average Lx/SFR ratio 
during the starburst phase is then a factor of ^ 2 lower 
than that of the normal SF population (in the common 
mass range, Fig.l lower panel). Keeping in mind that 
Lx and SFR are dominated by different populations, one 
possible interpretation of this result could be a delay in 
the BHAR enhancement relative to that of the SFR, as 
suggested by Hopkins et al. (2012, see also Lapi et al. 
2014). However, as discussed below, our observations are 
inconsistent with the predictions of this model, and, thus, 
a simpler argument is that mergers are more efficient at 
boosting SFR than BHAR. On the other hand, Alexan¬ 
der et al. (2008) found a similar lagging between the 
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Main statistical properties of the sample 


Main Sequence 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

log(Mass range) 

N stacked 

^detections 

< M, > 

< SFR > 

< L 2 - i 0 keV > 

10.09-10.42 

3215 

56 

10.24±0.09 

37.1±7.5 

''^-0.34 

10.42-10.75 

2158 

119 

10.56 ±0.09 

65.3±13.1 

O.UD_o.46 

10.75-11.25 

984 

124 

10.88±0.12 

149.3±29.9 

14 20'^°'®^ 

11.25-12.00 

36 

10 

11.31±0.11 

310.8±62.2 

cc 1 n + 10.96 

Starbursts above MS 

log(Mass range) 

-N stacked 

^detections 

< M, > 

< SFR > 

L2-10fceV 

10.25-10.75 

48 

12 

10.52 ± 0.15 

484.0 ± 235.4 

13.0-19.4 

10.75-11.25 

25 

6 

10.89 ± 0.13 

749.6 ± 328.8 

28.1-40.9 

Quiescent below MS 

log (Mass range) 

-N stacked 

^detections 

< M. > 

< SFR > 

L 2 - 10 keV 

10.42-10.75 

381 

4 

10.62 ± 0.09 

- 

0.63-2.88 

10.75-11.25 

627 

20 

10.95 ± 0.13 

- 

2.37-4.23 


TABLE 1 

Notes: (1) AL* range in solar units. (2) Number of X-ray stacked sources. (3) Number of sources detected in the Chandra 

X-RAY catalogues. (4) AVERAGE STELLAR MASS, ERRORS ARE THE 1-a OF THE MASS DISTRIBUTION IN EACH BIN (Ioq^Mq)). (5) AVERAGE 

SFR (MQ/yr). (6) Average intrinsic AGN rest-frame 2-10 keV X-ray luminosity (10'^‘^ergs~^cm.~^) obtained by combining 

STACKED AND DETECTED X-RAY SOURCES. WHERE THE SINGLE STACKED X-RAY FLUX IS NOT DETECTED, WE REPORT 2o- UPPER AND LOWER 

LIMITS. 


SFR and the BHAR in a sample of X-ray-obscured Sub¬ 
millimeter Galaxies at z ^ 2. Moreover, state of the art 
semi-analytic models, based on a merger-driven scheme 
for starbursts and AGN (Lamastra et al. 2013), cor¬ 
rectly foresee the observed increase of Lx at increasing 
distance from the MS (i.e. the starburstiness) but they 
predict the same enhancement of both Lx and SFR. 

Thanks to the improved statistics of our sample 
at higher masses, we found that the L^/SFR ratio 
is not invariant with the stellar mass, but scales as 
\og{Lx/SFR) oc MO'43 ±o o 9 (Fig. 1, lower panel). By 
contrast, M12 reported a constant value for this ratio as 
a function of mass. Since Lx traces the BHAR (Mbh), a 
constant Lx/SFR ratio translates into a Mbh/M^, con¬ 
stant ratio {Mbh oc M*, see formalism in M12). By inte¬ 
grating the superlinear relation reported here, we instead 
find that Mbh oc M^-^, implying a steepening of the 
Mbh/M^ distribution at z ~ 2 (in line with Kormendy 
& Ho 2013, who find a mass dependence of Mbh/M tuige 
in the local Universe). However, we note the Lx/SFR 
trend presented here is consistent to within the errors 
bars shown in M12 within the common mass range (see 
cyan points in the lower panel of Fig. 1), and it is 
also consistent within the observed scatter of the local 
Mbh/M^, relation (see Lasker et al. 2014). Our results 
also suggest that BHs accrete more mass with respect 
to their host galaxies at the high mass end, as also seen 
in a small sample of luminous, obscured QSOs at z ^ 2 
(Bongiorno et al. 2014). 

The definition of starbursts in our selection is consis¬ 
tent with what is envisaged in hydrodynamical simula¬ 
tions of merger galaxies that include feedback from a BH. 
This is shown in the panel a) of Figure 3, where we re¬ 
port the SF history of the Di Matteo et al. (2005) and 


Hopkins et al. (2012) models, normalized to the pre¬ 
merger phase. These models predict a sudden enhance¬ 
ment of Lx (panel b) and Lx/SFR (panel c) during the 
merger phase, up to a factor of 100. Within our sam¬ 
ple of starburst-, sSFR-selected galaxies we do not find 
an enhanced Lx/SFR ratio, which is actually slightly 
sub-MS. Clearly, the mere fact that these galaxies are 
starbursting indicates that the current AGN activity is 
not able to suppress SF in these (merging) systems. On 
the other hand, the peak of BHAR may be delayed with 
respect to the peak SFR (e.g., Chen et al. 2013, Hickox 
et al. 2014), taking place when the starburst has al¬ 
ready subsided, which is however contrary to the above 
simulations, where starburst and rapid BH growth are 
simultaneous events. 

We caution that among SF galaxies (for MS and in 
particular for the starbursts, Lanzuisi et al. 2014) we 
might have missed a fraction of Compton Thick AGN 
for which the adopted conversion factors to derive the 
intrinsic Lx might be underestimated. Surveys sensitive 
to more penetrating hard X-rays (e.g. from NuSTAR) 
are required to unravel this hidden extinction (Mullaney 
et al. in prep.). 
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